Thermodynamics with density and temperature dependent particle masses and 
properties of bulk strange quark matter and strangelets 
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I. INTRODUCTION 

One of the most exciting possibility from QCD is that 
hadronic matter undergoes a rich and varied phase land- 
scape with increasing densities. At extremely high den- 
sities, the mass of strange quarks becomes unimportant 
and all the three flavors of it, d, and s quarks can be 
treated on an equal footing. Consequently, quark mat- 
ter is in the color-flavor-locked (CFL) phase [l| and/or a 
new gapless CFL phase (gCFL) 2]. However, if densities 
are not that high, the strong interactions between quarks 
become important, and quark matter is in the unpaired 
phase. Presently, RH1C is teaching us about the prop- 
erties of the hot but not asymptotically hot quark gluon 
plasma 0,0. Actually the future FAIR project which 
will be built in the coming years at GSI in Germany is 
targeted towards the physics of ultradense matter like it 
is found in neutron stars [f| . 

The original idea of Witten is that strange quark mat- 
ter (SQM), rather than the normal nuclear matter, might 
be the true ground state of the strong interaction |g. 
Immediately after Witten's conjecture, Farhi and Jaffe 
showed that SQM is absolutely stable near the normal 
nuclear saturation density for a wide range of parame- 
ters 0- Now SQM has been investigated for more than 
two decades since the pioneer works of many authors 
0, 0, H, @ • Because the lattice gauge theory still has dif- 
ficulty in the consistent implementation of chemical po- 
tential in numerical simulations presently [T(| while the 
perturbative approach is unreliable at the strong coupling 
regime, phenomenological models reflecting the charac- 
teristics of QCD are widely used in the study of hadrons, 
and many of them have been successfully applied to in- 
vestigating the stability and properties of SQM. 

An important question in the study is how to incorpo- 
rate pressure balance. Basically, one uses a model to give 
the thermodynamic potential, then add to it a constant 
to get mechanical equilibrium. This is the famous 'bag' 



mechanism. Many important investigations have been 
done in this direction [lllIlHQEHIIlHIl! 

A common feature of these investigations is 
that quark masses are constant, so the normal thermody- 
namic formulas can be used without thermodynamic in- 
consistency problems. Actually, however, it is well known 
in nuclear physics that particle masses vary with environ- 
ment, i.e., the density and temperature. Such masses are 
usually called effective masses |23|, |2J, |25|, |2(j . Effective 
masses and effective bag constants for quark matter have 
been extensively discussed, e.g., within the Nambu-Jona- 
Lasinio model [27j and within a quasiparticle model [2^| . 
In principle, not only masses will change in the medium 
but also the coupling constant will run in the medium 
|29|. The question now is how to treat the thermody- 
namic formulas which do not violate the fundamental 
principles of thermodynamics, when introducing density 
and/or temperature dependent masses. In fact, a lot of 
problems have been caused in this field. 

There exist in literature several kinds of thermody- 
namic treatments. The first one uses all the thermody- 
namic formulas exactly the same as in the constant-mass 
case [3(3, 0, . The second treatment adds a new term, 
originated from density dependence of quark masses, to 
both pressure and energy |33| . These two treatments 
were later proved to be inconsistent with the necessary 
thermodynamic requirement: their free energy minimum 
do not correspond to zero pressure [34| . The third treat- 
ment adds the term from the density dependence of quark 
masses to the pressure but not to the energy, and the 
inconsistency disappear |34j . Another treatment is the 
addition of a new term to the thermodynamic potential. 
This has been done when masses depend on either tem- 
perature [35| or chemical potential/density [28|,|3(|. How- 
ever, if particle masses depend on both chemical/density 
and temperature, this way meets difficulties. We will 
discuss it further in the next section. 

With the third thermodynamic treatment and the cu- 
bic root scaling Zheng et al. have studied the viscos- 
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ity of SQM and calculated the damping time scale due to 
the coupling of the viscosity and r mode [33] • This model 
has also been applied to investigating the quark-diquark 
equation of state and compact star structure [39l |. 

Another important progress has been made recently by 
Zhang et al. 0, EH O E3| ■ They extended the quark 
mass density dependent model to finite temperature to 
let the model be able to describe phase transition. In 
this case, the quark masses depend on both density and 
temperature and the permanent confinement is removed. 
They use the third thermodynamic treatment mentioned 
above [34| and parametrize the common interacting part 
of quark masses as m i = ^o(3^b) _1 [l — (T/T c ) 2 ] 
with T c being the critical temperature. Later, they found 
this parametrization causes an unreasonable result: the 
radius of strangelets decreases with increasing temper- 
ature 41]. So they added a new linear term, and the 
parametrization became |4lt |42| 



mi = 



Bp 
3n B 



T , T 
a h b — 



(1) 



This extension of the model has soon been a ppli ed to 
the study of strangelets 0, E2> 0] > dibaryons |44| , and 
proto strange stars 0]i e * c -- 

The purpose of the present paper is two-folded. First, 
we would like to point out that the thermodynamic 
derivation in Ref. |34j is mainly concentrated on den- 
sity dependent masses. When masses are also tempera- 
ture dependent, there are special issues to be considered. 
We will prove that the temperature dependence of quark 
masses causes another term which contributes to the en- 
tropy and energy. With the new term, originated from 
the temperature dependence of masses, added to the en- 
tropy, the quark mass scaling Eq. Q has a serious prob- 
lem, i.e. limT-^o dmi/dT ^ which leads to a non-zero 
entropy at zero temperature, violating the third law of 
thermodynamics. Therefore, our second purpose is to de- 
rive or suggest a new quark mass scaling, and then study 
the properties of bulk SQM and strangelets. An inter- 
esting new observation is that low mass strangelets near 
(3 equilibrium are multiquark states with an anti-strange 
quark, such as the pentaquark (u 2 d 2 s) for baryon number 
1 and octaquark (u 4 d 3 s) for dibaryon etc. 

The paper is organized as follows. In Sec. [HI we derive, 
in detail, the thermodynamic formulas at finite tempera- 
ture with density and/or temperature dependent masses, 
and show why, in the present version of the quark mass 
density and temperature dependent model, one should 
require 



lim = 0. 

T^O dT 



(2) 



In the subsequent Sec. IIIII a new quark mass scaling at 
finite temperature is derived or suggested based on chi- 
ral and string model arguments, and accordingly in Sec. 
IIVI and Sec.0 respectively, the thermodynamic proper- 
ties of bulk SQM and strangelets at both zero and finite 



temperature are calculated with the new thermodynamic 
formulas and the new quark mass scaling. Finally, a sum- 
mary is given in Sec. I VII 



II. THERMODYNAMICS WITH DENSITY AND 
TEMPERATURE DEPENDENT MASSES 

Quasi particle models have been explored in great de- 
tail over the past 10-15 years |46|. It is well understood 
how to construct thermodynamically consistent models 
when the masses depend on the chemical potential and 
temperature. In the present model, however, there are 
three differences. First, the particle masses depend on 
density and temperature, not directly on chemical po- 
tential and temperature. Secondly, the density and tem- 
perature dependence is independently determined from 
other arguments. And thirdly, the finite size effects have 
to be included. Therefore, we derive, in the following, the 
thermodynamic formulas suitable for the present case. 

Suppose the thermodynamic potential is known as a 
function of the temperature T, volume V, chemical po- 
tentials {(J-i}, and particle masses {m;}, i.e., 



n = n(T,{ tH },{m i },v). 



(3) 



We here explicitly write out the arguments to make the 
meaning of partial derivatives clear in the following. If 
the masses m, are constant, other thermodynamic quan- 
tities can be obtained from normal formulas available in 
textbooks. Here the quark masses are density and/or 
temperature dependent, i.e., 



nii = mi(n h ,T), 



(4) 



where the baryon number density is connected to the 
particle numbers {iVj} and volume V by 



"b 



Ev 3 



with rii 



(5) 



To study this question, we start from the fundamental 
derivation equality of thermodynamics, i.e., 



dE = TdS - PdV + E VidNi. 



(6) 



This is nothing but the combination of the first and sec- 
ond laws of thermodynamics. It means that the energy 
E is the characteristic function, i.e., all other quanti- 
ties can be obtained from it, if one takes the entropy 5, 
the volume V, and the particle numbers {Ni} as the full 
independent state variables. But it is sometimes conve- 
nient to take (T, V, {Hi}) as the full independent state 
variables. In this case, the characteristic function is the 
thermodynamic potential Cl which is defined to be 



fl = E - TS - E f^Ni 
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because adding —d(TS + J2i HiN) to both side of Eq. 
© will give 



dVL = -SdT - PdV - Nid^i 



(8) 



Another important characteristic function is the free 
energy F. It is defined to be 

F = E-TS. (9) 
Then the corresponding basic derivation equation is 

dF = -SdT - PdV + Y VidNi (10) 



which can be obtained by adding —d(TS) to Eq. 
©. Therefore, the independent state variables are 
(T, V, {A}) in this case, i.e., 



F = F(T,V,{N i }). 



(11) 



According to the second term on the right hand side of 
Eq. I|1(J|) . one has a general expression for the pressure 



P 



dF 
dV 



T,{ N i} 



(12) 



Here F should be expressed as a function of (T, V, {N}), 
and the derivative is taken with respect to the volume at 
fixed T and {N}. Comparing Eqs. 10 and leads to 
the basic relation between thermodynamics and statis- 



F = ft + ^2 mNi 

i 

Substituting this into Eq. {T2J| gives 



(13) 



P : 



dn 

~d~V 



T,{Ni} 



av' 



(14) 



Because the independent state variables here are 
(T, V, {Ni}), the chemicals {&} in ft [see Eq. ©] should 
be expressed as a function of (T, V, {A/j}), i.e., 



m = m(T,v,{Nk}). 



(15) 



So the total derivative of ft with respect to V at fixed T 
and {Ni} is 



dft 
dV 7 



T,{N k } 



an 
av 



dm av + ^' 



dm; dV 



where 
dm; 



dm; 



av 



av 



dmi drib drrii 



T {JVfc} dn h dV V dn b 



(16) 



(17) 



Consequently, substitution of Eq. (|16fl into Eq. I|12|) gives 



„ 9ft rih v 9ft 9m, 



9^ y 



(9mi 9nb 



Similarly, for the entropy, we have 



9ft \ d/j,i 

dJUJ av" 
(18) 



S 



dF 
dT 



V,{N k } 

Substitution of 



dft 
dT 



V,{N k } 



E^> 



djM 
1 dT' 



dft 
dT 



9ft x -> 9ft 9/i.; x -> 9ft 9m,; 



9m 2 dT 



leads to 



9ft \ 9ft 9m, x / - 



dT ^ dm, dT 

i 

The energy can be obtained from Eq. J7J 
E = ft + ^ /^JVi + T5. 



9ft \ d/ii 
dJTj W 



(19) 



(20) 



(21) 



(22) 



Replacing S here with the expression in Eq. (|2~TJ) . one has 



- x - - m 9ft 9ft 9m,- 

£ = ^ + E ~ T Wr " T E ' 



9T 

9ft \ 9/Lti 
9^J 9T' 



9m, 9T 



(23) 



Now we need iVf (number of flavors) equations to 
connect T, V, {^i}, and Ni, so that the functions 
Hi(T, V, {Nk}) in Eq. (|15|l can be obtained. Presently, 
nearly all relevant models adopt [2?| |3(3, |U 0, HH HE 

USEES! El EH 



N 



9ft 
9/ij 



T,V,{m fc } 



(24) 



Please note, Ni also appears on the right hand side of 
this equation through m,k = mfc^ 7Vj/[3V], T). Then 
/Xi can be solved as a function of T, V, and {A^} from 
these equations. 

If Eq. applies, the last term in Eqs. fl%|) . an d 
(|23|l vanishes. So all the formulas will take the simplest 
form. If define ft = fl/V, E = E/V, S = S/V, m = 
Ni/V, and use V = (4/3)7ri? 3 , then Eqs. JTSJ, jgj) , 
lfT5)l. and Eq. become, respectively, 



^ ^ dmi 

dV dmi dnb 



(25) 



^ 9ft 9ft ^ 9ft 9m; 



S = 



9ft 

"df" 4- 



9T 

9ft 9mi 
9m, 9T ' 



(27) 



4 



F 



dn 



(28) 
(29) 



Compared with thermodynamic formulas in the nor- 
mal case, both the entropy and energy have a new term 
from the temperature dependence of the masses, while 
the pressure has a new term due to the density depen- 
dence of the masses. The second term in Eq. I|25f) exists 
when the finite size effect can not be ignored, no matter 
the masses are constant or not. These new terms are 
understandable with a view to the equalities 



dn 



V,n 



dn 



(30) 



7> 



and the total derivative rules in mathematics. Here one 
should pay special attention to the difference between the 
total derivatives and the partial derivatives in Eq. 118(1 
and Eq. (J3TJ. 

/From the viewpoint of quasiparticle models, the first 
and second terms in the pressure Eq. 1(2511 are the normal 
quasiparticle contributions, while the last extra term is 
the contribution of mean- field interactions. 

Normally, the first term on the right hand side of Eq. 
((27(1 goes to zero when the temperature approaches to 
zero, i.e., lim-r^o dn/dT = 0. However, the first factor 
of the second term does not, i.e., limx_,o dn/drrii ^ 0. 
Therefore, one must require 



,. dm l 
inn — — 

T-o dT 



= 



(31) 



to be consistent with the third law of thermodynamics, 
i.e., limj^o S = 0. Because Eq. 12711 depends on the 
model assumption Eq. Eq. l(3*Tl) is a model depen- 

dent requirement. Consequently, the thermodynamic for- 
mulas Eqs. ((25(l - i(29|) are merely valid for systems whose 
interactions meet this requirement. In the subsequent 
section, we will see that Eq. ((3111 is indeed fulfilled by 
a new quark mass scaling for the strong interactions of 
quarks. 

In Ref. [34f. we mentioned a necessary condition any 
consistent thermodynamic treatment must satisfy. At 
finite temperature, we also have a similar criterion: 



P=- 



dF 
dnu 



T,{N k } 



dV 



' dny 



(32) 



T,{N k } 



In obtaining the second equality of Eq. 1(32(1 . the equali- 
ties V = J2i Ni/{3n h ), dn h /dV =_-n h /V, and F = VF 
have been used. Because F/n^ — F / ATj/3) is the free 
energy per baryon, Eq. 1(32(1 shows explicitly that the free 
energy extreme occurs exactly at zero pressure. This is 
why people are interested in the free energy minimum, 
rather than the energy minimum, at finite temperature, 
to look for mechanically stable states. 

The extra term or the last term in Eq. ((25(1 is impor- 
tant even in the MIT bag model when one introduces a 



density dependent bag constant B(rib)- In this context, 
the extra term is n^dB / dn^. If it is not included, the zero 
pressure will not be located at the free energy minimum. 
In Ref. this term has been considered in the calcu- 
lation of hadron-quark phase transition in dense matter 
and neutron stars within the bag model. 

The derivation process of the pressure is also instruc- 
tive to the case when one wants to include the Coulomb 
contribution. In this case, the energy density, accord- 
ingly the thermodynamic potential density, gets a term 
Ecoui{{nk}, R)- Correspondingly the pressure gets 



Coul 



^Coul 



Coul 



RdE t 



Coul 



di 



3 dR 



(33) 



In literature, there is another approach to have ther- 
modynamic consistency by adding a term B* to the orig- 
inal thermodynamic potential density. This makes the 
total thermodynamic potential density becomes n + B* 
[28l l35l l3(| . We comment that this can be uncondition- 
ally done merely in two special cases, i.e., at finite tem- 
perature with zero chemical potential and at finite chem- 
ical potential with zero temperature. The expression of 
the added term for the former case is 1351 : 



B*{T) 



T on 

To dm 



7>=0 



dm 
~dT 



dT, 



while that for the later case is |2f 



^ dm 



dm 

-r-dfx. 



(34) 



(35) 



However, if one wants to extend this to both chemi- 
cal potential/density and temperature dependent masses, 
difficulties arise. To perform the integration in Eqs. 134() 
and 1(35(1. one should use m = m(T) to replace the m on 
the right hand side of Eq. 1(34(1 or apply m — m(/x) to 
the right hand side of Eq. l(3*5)l. If m, = mj(T, {/Ufc}), the 
directly extended integration is 



B*{T,n) = -J 



E 



dn dm, 



dT 



dn dm, 



dmj dT ' 4-^ dm* duj 

j J 

(36) 

In the above, (Tq,/j,q) is some reference point. Eq. 1(36|) 
is a multi-dimensional integration. To let it be path- 
independent, one must mathematically require Cauchy 
conditions. If all chemical potentials are equal, for ex- 
ample, the Cauchy condition is 



E 



d 2 n dmj 



d 2 n dmj 



drriidfx dT dmjdT dji 



0. (37) 



Such an example has recently been given in Ref. |48j . 
where the Cauchy condition Eq. 13711 . or the equivalent 
Maxwell relation Eq. (7) in Ref. 48], is fulfilled by solv- 
ing the equation for the coupling in the masses, As the 
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Eq. (8) of Ref. |4g indicated. However, if the masses 
are completely determined from other arguments, and so 
there are no adjustable parameters, this thermodynamic 
treatment may fail. When masses depend on density and 
temperature, rather than directly on chemical potential 
and temperature, the case becomes much more involved. 
Also, if finite size effects can not be ignored, or in other 
words, fl depends explicitly on the volume or radius, not 
merely the integration in Eq. (|36|l becomes much more 
difficult or impossible, the integration in Eq. (|34f> or in 
Eq. (|35[l has an unknown function of the volume or radius 
as well. 

These difficulties are not surprising. In fact, we have 
proved, from Eq. © to Eq. (23), that Eqs. igSJl-® are 
inevitable consequences of Eq. J3J with Eq. (|24[1 . In this 
paper, we will apply these formulas to the calculation of 
properties of both bulk SQM and strangelets. 



III. DERIVATION OF QUARK MASS SCALING 

In the preceding section, we have derived the ther- 
modynamic formulas suitable for systems with density 
and/or temperature dependent masses. In this section, 
we derive quark mass scaling by a similar method as in 
Ref. [13. 

Let's schematically write the QCD hamiltonian density 
for the three flavor case as 



H 



QCD 



E 

q—u,d,s 



m qa qq 



(38) 



where m q o (q=u,d,s) are the quark's current mass, is 
the kinetic term, Hi is the interacting part. 

Now we want to include interaction effects within an 
equivalent mass m q . For this purpose we define an hamil- 
tonian density of the form 



#cqv = H k + ^2 m q qq, 

q—u,d,s 



(39) 



where m q is our equivalent mass to be determined. We 
firstly divide it into two parts, i.e. 



m„ = m 



q o + mj. 



(40) 



The first part m q o (q — u, d, s) are the quark current 
masses while mi is a common part for all the three flavors 
to mimic the strong interaction. Obviously we must re- 
quire that the two hamiltonian densities H oqv and Hqqd 
have the same eigenenergy for any eigenstate |\E'), i.e. 



(*|ff cqv |*) = (*|Uqcd|*} 



(41) 



Applying this equality, respectively, to the state \nb,T) 
and the vacuum |0), and then taking the difference, we 
have 



(H c , 



(H 



cqv/Oj ' 



(H 



QCD/ rib ,T 



(ff QCD )o, (42) 



where the symbol definitions (0) nb ,r = (n>b, T\0\n\>, T) 
and (O)o = (OjOjO) have been used for an arbitrary op- 
erator O. Then solving for m\ from this equation gives 



mi 



(Hi) 



E Km)n b ,T - {qq)o\ ' 

q— u,d,s 



(43) 



where (Hi) = (Hj) nb! T — (Hj)q is the interacting part 
of the energy density from strong interactions between 
quarks. It can be linked to density rib and temperature 
Tby 



(Hi) =3n b v(f,T). 



Here 



"(- 



1/3 



(44) 



(45) 



is the average distance of quarks at the density n^, v(f, T) 
is the interaction between quarks at density rib and tem- 
perature T. Because we are interested in the confinement 
effect, while the lattice simulation ^] and string model 
investigation show that the confinement is linear, we 
write 



v(n b ,T) = a(T)r. 



(46) 



The temperature dependence of the string tension 
<j(T) can be obtained by combining the Eqs. (94) and 
(91) in Ref. 0: 



<j(T) = (T exp — 

a \ 1 



(47) 



where a is the lattice spacing while do is the string ten- 
sion at zero temperature. The value of cto from potential 
models varies in the range of (0.18, 0.22) GeV 2 |5^ . 

For convenience, let's define a dimensionless constant 
A = 2aoa/T c with T c being the critical temperature. 
Then substituting a = \T c /(2a ) into Eq. (S) gives 



cr(T) = (T 



f 



8T ( T c 



(48) 



Because the string tension should become zero at the 
deconfincment temperature, the value of A is determined 
by the equation ct(T c ) = whose solution is 



A = LambertW(8) w 1.60581199632. 
Accordingly, Eq. 146fl becomes 



v(n h ,T) = 



1/3 



8T 

^atT 6 



— AT C /T 



(49) 



(50) 



The inter-quark potential has been also studied by 
comparison to lattice data in Ref. [5i|- Replacing the 
factor exp(— jir) in the Eq. (2.8) there with I — fir, one 
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FIG. 1: Comparison of different temperature factors. The 
full line, dashed line, and dotted line are, respectively, for 
l-^e- ATc/T , 1- (£) 2 , and 1-0.65^-0.35 (£) 2 , where 
T c is the critical temperature and A = LambertW(8) ~ 1.6. 



will find the linear confining part. The temperature fac- 
tor used there is 1 — (T/T r ) 2 and this temperature factor 
has also been used in Ref. [4(|. However, because of a 
problem with the radius-temperature relation mentioned 
in the introduction part, the factor has been modified to 
1-0.65T/T C -0.35(T/T C ) 2 [U |H |H |3 . The tempera- 

8T 



ture factor derived in the present paper is 1 — -^e XT "/ T . 



These temperature factors are compared in Fig. Q 

Chiral condensates have been extensively studied in lit- 
erature [54l IHlj . In principle, they depend on both den- 
sity and temperature. Presently for simplicity, we only 
consider their density dependence and use the model- 
independent result [26ll56| 



with 



(qq)o p* ' 



(51) 



(52) 



When taking 140 MeV for the pion mass m n , 93.3 MeV 
for the pion decay constant /„-, and 45 MeV for pion- 
nucleon sigma term on, one has p* ~ 0.49 fm -3 . 

Substituting Eqs. piTJl. (JSUJl, and g3J to Eq. we 
obtain 



mi(n h ,T) 



D 



1 



8T 

at: 



exp 



-/V 



T 



(53) 



Here z = 1/3. Please note, many density and tempera- 
ture independent constants such as the vacuum conden- 



sates, the string tension at zero temperature, and p* et 
al. are grouped into a constant D, i.e., 



D 



3(2/7r) 1 / 3 <Jop* 



(54) 



Taking a = 0.18 GeV 2 and p* = 0.49 fm~ 3 , \f~D value 
is in the range of (147, 270) MeV when — J2 q (QQ)o varies 
from 3x (300 MeV) 3 to 3x (200 MeV) 3 . Becasue of many 
uncertainties within the relevant quantities and the fact 
that Eq. H51(l is not exactly valid, at least it ignores tem- 
perature dependence and higher orders in density, we do 
not try to use the relevant quantities to calculate the 
value for D from Eq. I|54|) . Instead, we treat D as a free 
parameter to be determined by stability arguments, i.e., 
it makes the energy per baryon E/n^ at zero tempera- 
ture is greater than 930 MeV for two flavor quark matter 
in order not to contradict standard nuclear physics, but 
less than 930 MeV for three flavor quark matter so that 
SQM can have a chance to be absolutely stable. Obvi- 
ously, the range determined by this method depends on 
the thermodynamic formulas and the values of quark cur- 
rent masses. In the present calculation, we use m U Q = 5 
MeV, m d0 = 10 MeV, and m s0 = 120 MeV for the cur- 
rent masses involved. These conditions constrain \f~D to 
a very narrow range of (154.8278, 156.1655) MeV, and 
we take D = (156 MeV) 2 in this paper. 

The model described in the above is a combination of a 
chiral model and a string model. There should be noth- 
ing really wrong with it. In fact, its zero-temperature 
form has been successfully applied to studying the prop- 
erties of SQM @E3j calculating the damping time scale 
of strange stars due to the coupling of the viscosity and 
r mode |38| . and investigating the quark di-quark equa- 
tion of state and compact star structure 39]. Further- 
more, various applications of the conventional quasi par- 
ticle models with chemical and temperature de pen dent 
masses have turned out to be very successful [46|, |48| . 
Therefore, we will apply the specific model presented here 
to the investigation of SQM in bulk and strangelets in the 
following. 



IV. PROPERTIES OF BULK STRANGE 
QUARK MATTER AT FINITE TEMPERATURE 

As usually done in this model, the quasiparticle con- 
tribution to the total thermodynamic potential density 
of SQM is written as 



SI = 2_^fli(T, pi, mi 



(55) 



where the summation index i goes over u, d, s quarks 
and electrons. Anti-particles are treated as a whole with 
particles, i.e., the contribution of the particle type i to 
the thermodynamic potential density is 



a = - 
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p 2 dp, (56) 



where m„ /ij, and T are, respectively, the particle masses, 
chemical potentials, temperature, and = {p 2 +m 2 ) 1 / 2 
is the dispersion relation. The particle number density 
corresponding to the particle type i is obtained by n% = 
-dQ/dm, giving 
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The energy density is E = J2i Fi(T, fH, rrii) with 
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dp 
(58) 



The free energy density f, the entropy density S, and 
the pressure P are, respectively, 



F = £ F i = £ ( fii + ^) ' 
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(59) 



(60) 



(61) 



dn\, drrii 

In the above, the partial derivatives relevant to f2j are 
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(62) 



and 
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1 + eKp-MO/T 
+ ln[l + e -( e ^+«)/ T ] 

K P + Mi)/T ^ (g3) 



1 + e (ti,P+m)/T 



The masses of electrons/positrons (0.511 MeV) and 
neutrinos / anti- neutrinos (if any) are extremely small. So 
they can be treated to be zero. For massless particles, 
the relevant integrations in the above can be carried out 
to give 



Qi ~ -24{^ + 2 ^ T + 15" T 



(64) 
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At zero temperature, we have the familiar results 
9 



(65) 
(66) 
(67) 
(68) 
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(69) 
(70) 

(71) 



m 2 In 



iMil + ^Mi - ™?i 



(72) 



The above formulas show that the number density is 
an odd function of the corresponding chemical potential, 
i.e., particle and anti-particle numbers are opposite in 
sign. But other quantities, such as the thermodynamic 
potential, entropy, energy, and free energy, are all even 
functions. 

Suppose weak equilibrium is always reached within 
SQM by the weak reactions such as 

d, s ^ u + e + i7 e , s + u <-> u + d. (73) 

Correspondingly, relevant chemical potentials satisfy 

Md = M«: (74) 
+ = M« + Me- (75) 

We also have the baryon number density equality 

- (n u + n d + n s ) = n b (76) 
and the charge neutrality condition 



2«« - 2«d - ^n s - n e = 0. 



(77) 



Neutrinos are assumed to enter and leave the system 
freely. So their chemical potential /i„ is zero. From Eqs. 
(|64|l - (|68|) . they have no contribution at zero temperature, 
but contribute at finite temperature. 
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FIG. 2: The entropy per baryon of SQM as an function of 
temperature. It increases with increasing temperature. When 
the temperature approaches to zero, the entropy goes to zero, 
no matter the density is high or low. 



FIG. 3: Temperature dependence of the energy and free en- 
ergy per baryon of SQM. The energy is an increasing func- 
tion of the temperature while the free energy decreases with 
increasing temperature. 



For a given T and rib, the quark masses are obtained 
from Eqs. (gO) and J53J): 
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8T 
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(78) 



The corresponding partial derivatives are easy to get 
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(79) 



The chemical potentials fj,i (i — u, d, s, e) are obtained by 
solving Eqs. H74 M 77 )) . au other thermodynamic quanti- 
ties can then be obtained. 

Figure [21 shows the entropy per baryon as a function 
of temperature for different densities. It is an increas- 
ing function of the temperature and goes to zero at zero 
temperature. This is ensured by the fact that we have 
liniT^o dm q /dT = from Eq. Q8U[). It is interesting to 
note that we did not require this in deriving the scal- 
ing Eq. (|75|l in Sec. IIIII We got this automatically and 
naturally. 

Figure |3| gives the temperature dependence of the en- 
ergy and free energy. It is obvious that the energy is an 
increasing function of temperature while the free energy 
decreases with increasing temperature. 

In Fig. 21 we plot the density dependence of the en- 
ergy and free energy per baryon at different tempera- 
ture. The free energy minimum corresponds exactly to 
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FIG. 4: Density dependence of the energy and free energy 
per baryon of SQM at different temperature. On every free 
energy lines, the zero pressure point is exactly located at the 
minimum. However, they are generally not at the same points 
for the energy per baryon, except zero temperature where the 
energy and free energy is identical. 
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FIG. 5: Temperature vs baryon number densities for differ- 
ent pressure values. All lines go to the deconfinement tem- 
perature T c as the density approaches to zero. But at zero 
temperature, the densities corresponding to different pressure 
are different from each other, lower pressure corresponding to 
lower density. 



the zero pressure, satisfying Eq. I|32|) . However, these 
two points (zero pressure and the minimum) are gen- 
erally not the same for the energy per baryon at finite 
temperature. But at zero temperature they coincide be- 
cause the energy and free energy are equal to each other 
at zero temperature. 

As pointed out in Ref . :4C| , a density and temperature 
dependent mass model should have the ability of investi- 
gating phase transition. This is demonstrated in Fig. 
The figure is plotted by adjusting, for a given density, the 
temperature to such a value that it gives a definite pres- 
sure indicted in the legend. It is obviously shown that 
all lines go to the deconfinement temperature T c as the 
density approaches to zero. But at zero temperature, the 
density is different for different pressure. Higher densities 
correspond to higher pressure. 



V. PROPERTIES OF STRANGELETS 

To study strangelets, the special problem is to include 
the finite size effect. We do this by applying the multi- 
reflection method, originally comprised by Balian and 
Bloch [57]], later developed by Madsen 0, Farhi, Berger, 
and Jaffe 0, El etc., and applied to the mass density 
and temperature dependent model by Zhang and Su [42, 
liM |44| . We express the quasiparticle contribution to the 
thermodynamic potential density of strangelets as fl = 



J2i Mi(T, Mi, mi, R) with 

<>, - T ! {in [l + e-<»^>/ T 



In 



+ e -(Vf 2 +™f+^)/T] J n / dp (81) 



where the density of state n^(p, m-i, R) is given in the 
multi-expansion approach |57| by 



n-(p,m i; i?) = g. t 



3p 



6 



^2 + ^/s fa) + ^2 fc 



(82) 



Here Xi = rrii/p, the functions fs{xi) 0E1 and fc(xi) 
El are 



fs(xi) = -^-arctan^) 



and 



fc(xi) = 



12tt 2 



1 arctan(xi) 

2Xi 



(83) 



(84) 



Then all other thermodynamic quantities are straightfor- 
ward from Eqs. (|25[l - i|29[l . Here is the number density for 
flavor i: 



" = / [Vi -Vi \ ni(p,mi,R) dp, (85) 
h 

where ijf is the fermion distribution function, i.e., 



1 + e {V^+^T^i)/1 



(86) 



Because we treat the particles and anti-particle as a 
whole, the number densities rii can be both positive and 
negative theoretically. A negative particle number means 
anti-particles. Fig.[S]shows the chemical potential depen- 
dence of the number density for various temperature and 
mass. In general, the thermodynamic potential density 
is an even function while the number density is an odd 
function of the chemical potential, and the number den- 
sity is zero at zero chemical potential. If finite size effects 
can be ignored, as in Eq. (joTjl of the preceding section, 
this function is monotonically increasing, and so positive 
chemical potentials correspond to positive number den- 
sity. When the finite size effects are included (the surface 
term an d curvature term 0), the case becomes 

more complex. When the temperature is very high, the 
function is still monotonic. However, when the tempera- 
ture becomes lower than some special value, the function 
becomes non-monotonic, and the number density is neg- 
ative for some special positive chemical potentials. For 
smaller masses, the chemical potentials are smaller. So 
light quarks do not fall into this region of chemical po- 
tentials in actual cases. However, when quark masses 
become bigger, the corresponding chemical potentials be- 
comes also bigger. It is therefore possible that quarks 
with comparatively bigger mass happen to have positive 
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FIG. 6: Chemical potential dependence of the particle num- 
ber density at different temperature and mass. The radius is 
fixed to be 1 fm. 



FIG. 7: Chemical potentials and energy per baryon from 
bulk strange quark matter to finite baryon number at zero 
temperature and zero pressure. 



chemical potential and negative number density in some 
special cases. We will see such special examples a lit- 
tle later. One may concern that non-monotone of the 
number density leads to dnjd[i < which violates the 
stability condition. However, our results are always lo- 
cated in the regime where the derivative is positive. The 
regime is marked with full lines in Fig. El 
The free energy density is 



f = (At + fHrii, ) , 

i 

the energy density E is 



(87) 
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the entropy density is 

S = J2 
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and the pressure P is 
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In the above, the relevant partial derivatives are 
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1 + mf/p 2 



(94) 

Quark masses and relevant derivatives are still given 
by Eqs. |751)-(lgUl>. 

In Ref. [41j, charge neutrality is also imposed for 
strangelets. This is convenient for checking whether the 
formulas are continuous from bulk SQM to finite baryon 
number. Fig. Ogives the energy per baryon and chemical 
potentials versus baryon number at zero temperature and 
zero pressure. The horizontal lines are the corresponding 
values for bulk SQM. It is very clear that all quantities 
approach to the corresponding bulk values with increas- 
ing baryon number, and finite size effects destabilize low 
baryon number strangelets. 

Now we treat strangelets in another different way. In- 
stead of imposing the charge neutrality Eq. I|77|l . we re- 
quire that the electron and positron number densities are 
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FIG. 8: Energy and free energy vary with the radius at dif- 
ferent temperature. The mechanically stable radius (where 
the pressure is zero) is marked with a small circle on each line 
while triangles indicate the minimum. 



FIG. 9: Temperature dependence of the mechanically stable 
radius, energy and free energy per baryon for the baryon num- 
ber A — 20. Both the radius and energy per baryon increase 
with temperature. 



zero, i.e., 

n e = (95) 

because the radius of strangelets is much smaller than the 
Compton wavelength of electrons and positrons. Elec- 
trons and positrons are not involved in the strong inter- 
action, or in other words, they are not confined, so the 
finite size terms in Eq. I|82l) vanish for them. From Eq. 
|(SSJ, we have 

^ = y( T2 + ^i)- ( 9 6) 

Therefore, zero n e means zero With a view to the 
chemical equilibrium Eqs. I|74|l and J7SJ), we naturally 
get 

= /x<i = Ms- (97) 

In fact, Eq. (|97(l is the condition to find out the config- 
uration, which has the lowest energy per baryon, from 
the all possible strangelets with a fixed baryon number 
[l7j . Due to Eq. I|97|) . only one chemical potential is left 
independent. And it can then be determined by solving 

V nu + nd + ns) = Jm^m (98) 

for a given baryon number A, temperature T, and radius 
R. 



Figure [S] shows the energy and free energy per baryon 
as a function of the radius for A = 20 at different temper- 
ature. The points marked with an open circle are the me- 
chanically stable radius where the pressure is zero. The 
minimum of each line is marked with a triangle. Again we 
see that these two points are always the same on the free 
energy line. But they are different on the energy line at 
finite temperature. However, they coincide at zero tem- 
perature because the energy and free energy are equal at 
zero temperature. 

For a given A and T, the mechanically stable radius is 
obtained by adjusting it so that the free energy is mini- 
mized, or, simply by solving the equation 

P = 0. (99) 

The temperature dependence of the stable radius for 
A = 20 is plotted in Fig. [5] with a solid line. It is obvi- 
ously an increasing function of temperature. The corre- 
sponding energy per baryon is also plotted in the same 
figure, labeled on the right axis with a dotted line. It 
is also an increasing function of temperature. However, 
the free energy per baryon (dashed line) decreases with 
increasing temperature. 

Because charge neutrality is not imposed, strangelets 
here are charged. The electric charge can be calculated 
by 




(100) 



In Fig. 1101 we plot the charge to baryon number ratio 
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FIG. 10: Charge to baryon number of strangenets at finite 
temperature. At lower baryon numbers, the ratio exceeds 
1/2. This means a negative number for strange quarks, or 
the appearance of an anti-strange quark. The right axis gives 
the radius at T = 0. 



FIG. 11: Strangeness and energy per baryon versus baryon 
number at different temperature. For small baryon numbers 
at very low temperature, the strangeness per baryon becomes 
negative. 



as a function of the baryon number at different tempera- 
ture. This figure shows that the charge to baryon number 
ratio is a decreasing function of the baryon number and 
temperature. 

A noticeable feature is that the charge to baryon num- 
ber ratio at lower baryon numbers (A < 5 at zero temper- 
ature) becomes greater than 1/2. This is very different 
from that in the bag model where the charge to baryon 
number ratio is very small. For normal nuclei, this ratio 
reaches its biggest value 1/2. So it seems difficult to un- 
derstand a heavy positive charge at first sight. In fact, 
it is caused by the fact that f s is negative for 1 < A < 5 
and T = 0, i.e., anti-strange quarks, rather than strange 
quarks, appear. One can see this phenomenon clearly in 
Fig. 1111 where the free energy per baryon and the chem- 
ical potential have also been shown. Because of finite 
size effects, stranelets with very low baryon numbers are 
metastable for the parameters chosen. If one choose a 
bigger value for the current quark mass, the value for 
D should be smaller, and consequently, the mass for 
strangelets would be smaller. 

Fig. 1121 shows the quark configuration for low baryon 
numbers. To see the results clearly, we give the corre- 
sponding data in Tab. [I] The first column is the baryon 
number, the second to fourth column are, respectively, 
the quark numbers N u , Nd, and N s . Because shell ef- 
fects are not taken into account and beta-equilibrium 
is imposed, fractional quark numbers appear in Tab. 
1. Actual quark numbers should naturally be integers. 
So we approximate these real numbers to integers by 



int(ATj) + JVj/|JVj| (i = u, d, s; int means the number be- 
fore the decimal point). The results are shown in the 
fifth column. For A = 1, we have the pentaquark state 
(u 2 d 2 s _1 ). For A = 2, we have the dibaryon (u 4 d 3 s -1 ) 
or octaquark state. For A = 3, 4, and 5, we respec- 
tively have the multi-quark states (it 5 ci 5 ,s _1 ), (u 7 (i 6 s _1 ), 
and (it 8 d 8 s _1 ). A common feature of these states is that 
they all include an anti-strange quark. So we use 'slet' 
as the title of the fifth column. The charge number of 
these slets are given in the sixth column, while the sev- 
enth column gives the energies calculated by the present 
parameters (L> 1/2 = 156 MeV and m s0 = 120 MeV) with 
perfect (3 equilibrium. If one would like to produce 1540 
MeV (the actual + resonant mass) for A = 1, then one 
has to take D 1 ' 2 = 186 MeV and get 2856 MeV for A = 2. 
So we expect that the mass of the octaquark (u 4 d 3 s _1 ), 
if truely exists, is near 2856 MeV. For D 1 / 2 = 186 MeV 
and A > 3, the strange quark number becomes positive. 
So in this case we have only the pentaquark and the oc- 
taquark. Because of uncertainties in parameters, and also 
many other factors e.g. the perturbative interaction has 
not been included (the quark mass scaling is derived by 
assuming that the linear confinement interaction domi- 
nates), the concrete values should not be taken seriously, 
and further studies are needed. 

Recently, the pentaquark state + (154O) has aroused a 
lot of interest [58|. The width of 6+£l540) is very narrow 
with upper limit as small as 9 MeV [5!j • Cahn and trilling 
have extracted r(0 + ) = 0.9±0.3 MeV from an analysis of 
Xenon bubble chamber |60| . Although other hadrons like 
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FIG. 12: Quark configuration for low baryon numbers at 
T = 0. When 1 < A < 5, the s quark numbers are negative, 
indicating the appearance of anti-strange quarks. The stars 
are possible slets in the (A,S) plane. 



A 


N u 


N d 


N 3 


slet 


Z 


E (MeV) 


E* (MeV) 


1 


1.6951 


1.5416 


-0.2367 




1 


1290 


1540 


2 


3.2390 


2.9992 


-0.2382 




2 


2397 


2856 


3 


4.7570 


4.4360 


-0.1930 




2 


3471 


4131 


4 


6.2634 


5.8639 


-0.1273 


u'd b s-' 


3 


4528 


5385 


5 


7.7624 


7.2865 


-0.0489 


u^s' 1 


3 


5573 


6624 



TABLE I: Possible multi-quark states with an anti-strange 
quark. N u , Nd, and N 3 are the numbers of u, d, and s quarks, 
satisfying the weak equilibrium for the baryon numbers in 
the first column. The fourth column gives the possible slets, 
obtained by upgrading the N u , Nd, and N s to integers. The 
sixth column presents the electric charge number, while the 
seventh column is the corresponding masses with perfect j3 
equilibrium and D 1 ^ 2 = 156 MeV. The last column is for 
D 1/2 = 186 MeV, in which case the s quark number density 
for A > 3 becomes positive, namely, we have only the first 
two slets: the pentaquark and the octaquark. 

0(1020) and A(1520) are also narrow [T(0) = 4.26 ± 0.05 
MeV, T(A) = 15.6±1.0 MeV], we know what makes them 
narrow. However, we do not know, until now, why the 
+ should be so stable. We see from the above data that 
(u 2 d 2 s~ 1 ) is the multi-quark state which mostly satisfy 
the weak equilibrium for baryon number 1. This might 
serve as an explanation for the stability of 8 + (1540). 
Other slets, e.g. the dibaryon (w 4 (i 3 s _1 ), an octaquark 
state, should also exist, yet to be searched for by experi- 
ments though. 

In Ref. [15| , strangelets were also calculated to be heav- 
ily charged. But the electric charge is negative. There the 
investigation was concerned with how small metastable 



strangelets look like and might decay for different life- 
time. With the similar ideas, Ref. [i^] studied strangelets 
within density-and-temperature dependent quark masses 
and extending the findings in Ref. | l5| to finite tempera- 
ture. Present investigation concentrate on deriving ther- 
modynamic formulas and quark mass scaling, and find- 
ing the lowest-energy configuration from the strangelets 
with a fixed baryon number. Naturally, the observation 
of heavily positively charged strangelets, or multi-quark 
states with an anti-strange quark, depends on the value 
of the parameter D. If we took a much larger D value, 
the charge to baryon number ratio would also be small, 
or the anti-strange quark would not appear. However, 
bulk SQM has no chance to be absolutely stable in that 
case. 



VI. SUMMARY 

When masses are density and/or temperature depen- 
dent, the thermodynamical formulas are different from 
that for constant masses. We have derived a new set 
of thermodynamical formulas which can be used to cal- 
culate the properties of quark matter within a density 
and/or temperature dependent quark mass model. The 
new formulas are also instructive when one introduces a 
density and/or temperature dependent bag constant in 
the bag model etc. 

We have also argued for a new quark mass scaling at fi- 
nite temperature. The basic feature is that quark masses 
and their partial derivative with respect to the tempera- 
ture go to zero when the temperature approaches to zero. 
This ensures that all quantities restore to the density de- 
pendent model at zero temperature. It is especially im- 
portant that the entropy goes naturally to zero when the 
temperature approaches to zero, satisfying the third law 
of thermodynamics. 

With the new thermodynamical formulas and new 
quark mass scaling, we have studied the properties of 
bulk SQM and strangelets. It is shown that the free en- 
ergy minimum corresponds exactly to the zero pressure, 
both at zero and finite temperature. The mechanically 
stable strangelet radius increases with temperature. An 
interesting new observation is that low mass strangelets 
are heavily positively charged, or appear as multi-quark 
states with an anti-strange quark, such as the pentaquark 
(u 2 d 2 s) and the octaquark (u 4 d 3 s) etc., if bulk SQM is 
absolutely stable. 
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